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The magnetic field distribution around the vortices in TmN^EbC in the paramagnetic phase was 
studied experimentally as well as theoretically. The vortex form factor, measured by small-angle 
neutron scattering, is found to be field independent up to 0.6-ff C 2 followed by a sharp decrease at 
higher fields. The data are fitted well by solutions to the Eilenberger equations when paramagnetic 
effects due to the exchange interaction with the localized 4/ Tm moments are included. The induced 
paramagnetic moments around the vortex cores act to maintain the field contrast probed by the 
form factor. 

PACS numbers: 74.25.Op, 74.25.Ha, 74.70.Dd, 61.12.Ex 



The interplay between superconductivity and local 
magnetic moments is a fascinating problem, with rel- 
evance to a number of important unresolved questions 
such as the detailed nature of both high-T c and heavy- 
fermion superconductivity. Adding to this the usually an- 
tagonistic nature of superconductivity and magnetism, it 
is no surprise that materials which exhibit a coexistence 
of these two goundstates attract a lot of attention. 

The antifcrromagnctic members of the intermetallic 
nickclborocarbidc superconductors i?Ni2B2C (R = Ho, 
Er or Tm) have proved especially rich vehicles for such 
studies, displaying e.g. intertwined magnetic and super- 
conducting transitions as well as subtle changes in the 
superconducting characteristic length scales associated 
with the onset of antiferromagnetic ordering [H,0,0| ■ The 
exchange interaction TL s f = —I(gj — 1)J • s between the 
4/ localized moment J and the conduction electron spin 
s (gj is Lande g-factor and / is exchange integral) is im- 
portant in understanding systematic changes of both su- 
perconducting and magnetic transition temperatures Q . 
However, even in the paramagnetic state above the anti- 
ferromagnetic ordering temperature, Tjq, the conduction 
electron moment /x e = g^B s is subjected to an exchange 
field H cx = I(g.j — 1)3 /gus due to the field induced 4/- 
moments, yielding a "Zeeman" term Ti. s f = — H cx • \x e in 
the conduction electron Hamiltonian. 

Here we report on combined experimental and theoret- 
ical studies of TmNi2B2C, investigating specifically how 
the magnetic field profile around the vortices is influenced 
by the paramagnetic state. Using small-angle neutron 
scattering (SANS) we imaged the vortex lattice (VL) at 
several temperatures > Tn, and measured the magnetic 
field dependence of the form factor which reflects the field 
distribution around the vortices. In contrast to the usual 
exponential decrease with increasing field, the VL form 



factor in TmNi2B2C remains constant up to H ~ 0.6H C 2, 
followed by a sudden decrease as the upper critical field is 
approached. It is the striking departure from exponential 
behavior which is the central result of this Letter. 

The experimental results are compared to solutions of 
the quasi-classical Eilenberger equations, focusing on how 
the internal field distribution in the mixed state is af- 
fected by changes to the electronic vortex core structure 
due to the paramagnetism. Since M (H) is roughly linear 
below H C 2 [B[ the induced moment J is proportional to 
the applied field and can thus be treated as an effective 
Pauli paramagnetic effect, giving rise to a Zeeman en- 
ergy [iB where the parameter (i signifies the strength of 
the paramagnetic effect. The calculations show how the 
induced moments in and around the vortices grow with 
increasing applied field and thereby maintaining a high 
field modulation (and hence form factor), before they 
eventually spread out from the core region at high fields. 
The results of the calculations provide an excellent quan- 
titative agreement with the measured form factor. 

TmNi2B2C has a superconducting critical tempera- 
ture, T c = 11 K, and the Tm moments order antiferro- 
magnetically in a long-period transverse-modulated state 
below T N = 1.5 K [J, [1, @, 0|. At low temperature the 
magnetic moments are along the c axis, which conse- 
quently is the direction of the maximum magnetic sus- 
ceptibility Q . For magnetic fields applied parallel to 
the c axis H C 2 shows a non-monotonic behavior, reach- 
ing a maximum near T = 5 K and /io-ff C 2 = 1 T due 
to the Tm sublattice magnetization, decreasing upon ap- 
proaching the magnetic ordering temperature and finally 
increasing again below Tn 0, [H] • Previous studies showed 
simultaneous magnetic and VL symmetry transitions be- 
low Tn, as well as peaks in the FLL reflectivity associated 
with the magnetic transitions [l| . 
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FIG. 1: (Color online) Measured magnetic induction ver- 
sus applied field at 1.6 K for TmlS^EbC (squares) and non- 
magnetic LUM2B2C (circles). The top left inset show a VL 
diffraction pattern obtained at 0.2 T and 1.6 K. The bottom 
right inset show a schematic of the diffraction pattern, indi- 
cating the VL scattering vectors and opening angle. Open 
and closed circles represent peaks belonging to different do- 
main orientations, while x'es denote higher order reflections. 



SANS experiments were carried out at the Dll instru- 
ment at the Institut Lauc-Langcvin. The TmNi2B2C sin- 
gle crystal used in the experiment was grown using a high 
temperature flux method, using isotopically enriched n B 
to reduce neutron absorption [8[ . Incident neutrons with 
wavelengths of A rl = 6 - 8 A and a wavelength spread 
of AA„/A„ = 10% were used. The VL diffraction pat- 
tern was collected by a position sensitive detector. For 
all measurements, the sample was cooled in a horizontal 
magnetic field applied parallel to the crystalline c axis 
and the incoming neutrons. Measurements obtained at 
zero field were used for background subtraction. 

The VL was imaged as a function of field at temper- 
atures, T = 1.6 K, 3.5 K and 5.0 K. At all fields and 
temperatures a rhombic VL was observed, as shown in 
the insets to Fig. [T] The opening angle, (3, was found 
to decrease with increasing field indicating a continuous 
transition from a distorted square to a distorted hexag- 
onal symmetry in agreement with previous reports [l|. 
As a consequence of having a non-square VL pinned to 
an underlying square crystalline lattice, two VL domains 
were observed at all measured fields and temperatures. 

A direct measure of the magnetic induction, B, in the 
sample can be obtained from the VL scattering vectors. 
Using two scattering vectors belonging to the same do- 
main, the induction is given by 
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(1) 



where <fi = 20.7 x 10 4 TA 2 is the flux quantum. Fig. Q] 
shows the measured induction as a function of applied 
field for TmNi2B2C. To rule out the possibility of sys- 
tematic errors on the determination of B, mcasure- 
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FIG. 2: (Color online) Field dependence of the measured VL 
form factor in TmNi2B2C at 1.6 K. The lines show the form 
factor calculated using eqn. © with £ = 210 A and A = 780 
A (dotted line) or 600 A (dashed line). 



ments on non-magnetic LUM2B2C were performed im- 
mediately prior to the measurements on TmNi2B2C us- 
ing the same instrumental configuration. The measure- 
ments on LuNi 2 B 2 C yielded dB/d(^ H) = 1.003 ±0.006 
(dashed line) as expected for a non-magnetic supercon- 
ductor when H > H cl . For TmNi 2 B 2 C wc find B > fx H 
for the entire measured field range as seen in Fig. [l] in- 
dicating a significant paramagnetic contribution to the 
induction. Below - 0.6 T, dB/d((J, H) = 1.152 ± 0.004 
as indicated by the straight line in Fig. [TJ Taking demag- 
netization effects into account, this is in excellent agree- 
ment with magnetization measurements At higher 
fields B approaches [IqH, with the two fields seemingly 
merging at /zo-ff C 2 ~ 0.75 T. We do presently not have an 
explanation for the high field behavior of B. 

We now turn to the main focus of this Letter: Measure- 
ments of the TmNi2B2C VL form factor, F(q), which is 
the Fourier transform of the magnetic field modulation 
due to the vortices. Experimentally the form factor is 
related to the VL reflectivity by 



(2) 



where 7 = 1.91 is the neutron gyromagnctic ratio, t is the 
sample thickness, and q is the magnitude of the scattering 
vector 0| . Fig.[2]shows the VL form factor for TmNi 2 B 2 C 
at 1.6 K just above T^, obtained from the integrated 
intensity of the Bragg peaks, as the sample is rotated 
through the diffraction condition. 

Using the model obtained by Clem and valid for kql = 

A/£> 1 M- 



F(q) = B 



1 + AV ' 



g=V2a.q 2 + \- 2 ) 1/2 , (3) 



where K\ is a modified Bcsscl function, we have calcu- 
lated expected vortex form factors shown by in Fig. [5] 
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The dotted line corresponds to a penetration depth, A = 
780 A from literature [5j and a coherence length based on 
the upper critical field at 1.6 K, £ c2 = ^4>o/2ttH C 2 = 210 
A. As it is evident, the calculated form factor falls sub- 
stantially below our inferred values. Using a somewhat 
smaller value of the penetration depth, A = 600 A(dashed 
line), yields a better fit to the endpoints of the measured 
form factor, but fails to describe the non-exponential field 
dependence at all intermediate fields. Qualitatively sim- 
ilar results were also obtained at temperatures further 
above Tn as shown in Fig. [3] As we will show below 
the unusual field dependence of the form factor can be 
explained by a microscopic calculation taking into ac- 
count paramagnetic effects which modify the magnetic 
field profile around the vortices. 

The form factor F(q h k ) is calculated from the internal 
field distribution B(r) = J2 h , k F (<lh,k) exp(iq M ■ r) with 
the wave vector q h k = + fcq 2 , q x = (2n/a, —Tr/a y , 0) 
and q 2 = (2ir/a, ir/a y , 0), corresponding to VL unit vec- 



tors (a/2,a y ) and (a/2,- 



The intensity of the main 



peak at (h, k) = (1,0) gives the fundamental component 
F(q 1 ). To determine B(r) we selfconsistently calcu- 
late the spatial structure of the pair potential A(r) and 
the vector potential A(r) using the quasiclassical Eilen- 
berger theory in the clean limit [H [H, EE 0, El, El, 
including the paramagnetic contribution due to the ef- 
fective Zeeman effect through the exchange coupling of 
the conduction electron and TmNi2B2C sublattice mo- 
ments. The quasi-classical Green's functions g(u> n , k, r), 
/(w„,k, r), and /t(aj„,k, r) are calculated in the vortex 
lattice state by the Eilenberger equations 

{w n +kxB + k- (V + iA)}/ = A<f>(k)g, 
{co n + i/iB-k.(V-iA)}/ t = A*0*(k) 5 , (4) 

with g = (1 - // T ) 1//2 , Rc{g} > 0, the pairing func- 
tion <^(k), Matsubara frequency uj n = (2n + 1)ttT, and 
effective Zeeman energy /iB [16j. Here /i determines 
the strength of the paramagnetic effect. A simple two- 
dimensional Fermi surface is used, with a Fermi momen- 
tum unit vector given by k = (cos 9, sin 9) and < 9 < 
2ir. With the magnetic field applied along the z axis di- 
rection, the vector potential A(r) = x r + a(r) in 
the symmetric gauge, where B = (0, 0, B) is the average, 
uniform flux density and a(r) is related to the modulated 
internal field such that B(r)=B + Vx a ( r )- 

The selfconsistcnt conditions for A(r) and A(r) are 
given by respectively 



and 



A = . 9o A T J2 (r(k)(/ + / t> 

0<w;<w ctl t 



2T 

V x (V x A) = V x M para - — (klm{ ff }) k , (6) 



(5) 
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FIG. 3: (Color online) Comparison of measured and calcu- 
lated VL form factors in TmN^EbC at T = 1.6, 3.5, and 5.5 
K. The curves were calculated using the model described in 
the text, for T = 0.16T C and fi = 1.71 (A), T = 0.35T C and 
H = 1.28 (B), and T = 0.50T C and p = 0.86 (C). 



where 

k = v /7C(3)/8k G l 



k indicates the Fermi surface average, and 
kgl where kql is the Ginzburg- 
Landau (GL) parameter [ljj. In Eq. |JSJ), (g N a )~ 1 = 



lnT + 2T So<w i <w cut u i ' and we usc ^cut = 20fc B T c . In 
Eq. ([6]) both the diamagnctic contribution of supercur- 
rent in the last term and the contribution of the param- 
agnetic moment M para = (0, 0, M para (r)) with 



M para (r) 



pi 



are treated fully self-consistently [16j . As mentioned ear- 
lier, the vortices in TmNi2B2C at low temperature and 
intermediate fields form a distorted square VL [H, indi- 
cating a large fourfold anisotropy of the Fermi surface 
and pairing function (l8| . We consequently use a pairing 
function </>(k) = |v / 2cos26'| and a square VL configura- 
tion (a y = a/2). However, the overall qualitative features 
of the form factor do not depend much on these choices. 

The field dependence of the calculated |F(q = q 01 )| is 
shown in Fig. [3l where we have used HqH C 2(T = 0.5T C ) = 
1 T for comparison with the experimental data. The 
magnitude of |_F((7)| 2 depends on the GL parameter, and 
its gradient as a function of field is related to the para- 
magnetic parameter [i. At T — 0.5T C (5 K) values for 
K = 6.2 and fj, = 0.86 were chosen to obtain agreement 
between the calculated and measured form factor at low 
and intermediate fields. The small deviation close to H C 2 
may be due to the increasing deformation of the VL away 
from a square symmetry [l[ , which is not included in this 
calculation. While K is kept constant for the remainder of 
the calculations, the value of \i is expected to be propor- 
tional to the magnetization which in TmNi2B2C is dom- 
inated by the contribution from the Tm 4/- moments [H . 
The decrease of H C 2 below 5 K can thus be attributed 
to an increasing paramagnetic depairing. We therefore 
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FIG. 4: (Color online) Spatial structure of internal field B(r) within a unit cell of the square vortex lattice. Here B — 0.1 = 
0.36#c2(0.5T c ), ft = 6.2, T = 0.16T C , and = 0.02 (a) and 1.71 (b). A profile of the field distribution is shown in (c). 



determine values of [i such that they reproduce the sup- 
pression of H c2 (T), yielding (M = 1.28 at T = 0.35T C (3.5 
K) and /x = 1.71 at T = 0.16T C (1.6 K), corresponding 
to respectively H c2 (T)/H c2 (0.5T c ) ~ 0.85 and 0.75. As a 
consequence of the increasing value of the slope of the 
form factor at low fields changes from negative to positive 
(curves C — > B —> A in Fig. [3]) . As evident from Fig. [3] 
the calculated form factor at 0.35T C provides a good fit 
to the experimental data. At 0.16T C the calculated form 
factor captures the qualitative field dependence, but falls 
below the datapoints at low fields. The reason for this 
quantitative deviation is not clear, but it may be related 
to critical behavior due to the close proximity to Tn- 

The paramagnetic moments induced around the vortex 
cores enhance B(r), and consequently also the form fac- 
tor |_F(q 01 )| which acquires a paramagnetic component 
proportional to H . To visualize the contribution from the 
paramagnetic moments around the vortex cores, Fig. 0] 
shows the spatial structure of the internal field B(r) in 
a VL unit cell for /i = 0.02 (a) and 1.71 (b). The vor- 
tex field profile along the nearest neighbor direction is 
plotted in Fig. [UJc). This shows how the paramagnetic 
component is confined at the vortex center resulting in 
the enhancement of the internal field. 

Before concluding we would like to emphasize the rela- 
tive "simplicity" of TmNi2B2C as well as the theoretical 
model used here to describe the SANS results. In con- 
trast to the dr or (triplet) p-wave pairing observed in 
respectively the high T c 's and SrRuC>4 [19( the Cooper 
pairs in TmNi ^Bo C are singlets with only a modest 
gap anisotropy [20(. Likewise Pauli paramagnetic limit- 
ing and a possible non-uniform superconducting (FFLO) 
state which have recently received considerable attention 
in the heavy fermion superconductor CeCoIns [2lJ is not 
relevant in the case of TmNi2B2C 
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and maintaining the field distribution contrast of the VL. 
While our calculation used a simple model to describe the 
H- and T-dependences, it was still able to capture the 
qualitative and quantitative behavior of the form factor, 
emphasizing that paramagnetic effects are important in 
undcrstaning the vortex state in TmNi2B2C. 

After submission we became aware of a theoretical pa- 
per by J. Jensen and P. Hedegard (to be published in 
Phys. Rev. B) which also treats the anomalous field 
dependece of the form factor in TmNi2B2C. 

We are grateful to V. G. Kogan for numerous stim- 
ulating discussions. MRE acknowledges support by the 
Alfred P. Sloan Foundation. Work at Ames Laboratory 
is supported by the U. S. Department of Energy, Basic 
Energy Sciences under Contract No. W-7405-Eng-82. 



Instead we argue 

that one can consider TmNi2B2C as a "standard" para- 
magnetic (above Tn) superconductor, thus providing a 
very valuable reference for more exotic materials. 

In summary, we have presented combined experimen- 
tal and theoretical studies of vortices and the vortex 
lattice in TmNi2B2C in the paramagnetic phase above 
Tn- The physical picture which emerges is that the con- 
duction electron paramagnetic moments induced by the 
exchange interaction accumulate exclusively around the 
vortex cores, creating nano-tubes of Tm magnetization [22] 
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The orbital critical field in TmNi2B2C estimated from 
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{dH c2 /dT)\ Tc yields H£ h = 2.1 T 0, an order of meV HHp. 

magnitude smaller than the Pauli limiting field H^ 2 = 
V2A(0) /giiB ~ 17.7 T obtained using A(0) = 1.45 



